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We revisit the dyon-anti-dyon liquid model for the Yang-Mills confining vacuum discussed by Di- 
akonov and Petrov, by retaining the effects of the classical interactions mediated by the streamline 
between the dyons and anti-dyons. In the SU(2) case the model describes a 4-component strongly 
interacting Coulomb liquid in the center symmetric phase. We show that in the linearized screening 
approximation the streamline interactions yield Debye-Huckel type corrections to the bulk param¬ 
eters such as the pressure and densities, but do not alter significantly the large distance behavior 
of the correlation functions in leading order. The static scalar and charged structure factors are 
consistent with a plasma of a dyon-anti-dyon liquid with a Coulomb parameter P« 1 in the dyon- 
anti-dyon channel. Heavy quarks are still linearly confined and the large spatial Wilson loops still 
exhibit area laws in leading order. The t' Hooft loop is shown to be 1 modulo Coulomb corrections. 

PACS numbers: ll.15.Kc, ll.30.Rd, 12.38.Lg 


I. INTRODUCTION 

At asymptotically high temperature T, QCD-like the¬ 
ories are in a weakly coupled state known as the Quark- 
Gluon Plasma (QGP). In it semi-classical solitons - in- 
stantons and their constituents, monopoles etc - have 
large action S = 0{l/as) » 1. Their semi-classical 
treatment is parametrically reliable, but their density is 
exponentially suppressed by e“'®. Asa result their effects 
are small. 

However, as the temperature decreases the semi- 
classical action S decreases. Since the soliton density 
grows as a power of 1/T their contribution to the QCD 
partition increases. At a critical density fixed by Tc, con¬ 
finement sets in, and the near-zero expectation value of 
the Polyakov line (L) « 0 switches off the quark compo¬ 
nent of the QGP, as well as the (non-diagonal) gluons. 
Below the critical temperature Tc, the solitons dominate 
the field ensemble. 

The major questions at the transition point are: (i) Are 
these objects still made of strong enough fields, allowing 
for a semi-classical analysis; (ii) Are their interactions 
weak enough to preserve their individual identity; (iii) 
Are the semi-classical interactions in the thermal ensem¬ 
ble amenable to known methods of many-body theory. 
As we will argue below, two first questions will be an¬ 
swered in the affirmative, and the third also, provided 
the ensemble is dense enough. 

The instanton liquid model developed in the 1980's is 
an example of such a semi-classical treatment. In vac¬ 
uum at T = 0, the action per typical SU(3) instan¬ 
ton was found to be large with S ~ 12, and the inter- 
instanton and anti-instant on interactions tractable. The 
non-perturbative vacuum topological fluctuations are re- 
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lated to the explicit violation of the axial U(l), and the 
formation of fermionic zero modes. The collectivization 
of the fermionic zero modes leads to the spontaneous 
breaking of flavor chiral symmetry [1] (and references 
therein). More recently, instanton-induced effects were 
found to be important for hadronic spin physics [5]. 

However, around the critical temperature T ^ T^, in- 
stantons should know about the non-vanishing of the 
Polyakov line expectation value, also referred to as a non¬ 
trivial holonomy. Instantons with non-trivial holonomies 
were found in [3]. The key discovery was that large 
holonomies split instantons into constituents, the self¬ 
dual instanton-dyons. Since these objects have nonzero 
Euclidean electric and magnetic charges and source 
Abelian (diagonal) massless gluons, the corresponding 
ensemble is an “instanton dyon plasma” 

We remark that the electric dyon field is real in Eu¬ 
clidean space-time but imaginary in Minkowski space. 
The instanton-dyons are also referred to as instanton- 
monopoles or instanton-quarks. However, the notion of 
a non-zero holonomy and all the instanton-related con¬ 
tractions do not exist outside of the Euclidean finite-T 
formulation. On the lattice, both the electric and mag¬ 
netic charges of the instanton-dyons are observable by 
standard Gaussian surface integrals. 

Diakonov and Petrov |3] emphasized that, unlike the 
(topologically protected) instantons, the dyons interact 
directly with the holonomy field. They further suggested 
that since such dyon (anti-dyon) fields become signifi¬ 
cant at low temperature, they may be at the origin of 
a vanishing of the mean Polyakov line, or confinement. 
This mechanism is similar to the Berezinsky-Kosterlitz- 
Thouless-like transition of instantons into fractional “in¬ 
stanton quarks” suggested earlier by Zhitnitsky and oth¬ 
ers [6] , inspired by the fractionalization of the topological 
charge in 2-dimensional CPN models [7], although it is 
substantially different in details. It is also different from 
the random dyon-anti-dyon ensemble suggested earlier 
by Simonov and others [8]. It is not yet clear how this 
Euclidean mechanism relates to the the quantum con- 
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densation of magnetic monopoles suggested initially by 
t' Hooft [9] and Mandelstam [10], and subsequently sup¬ 
ported in the supersymmetric model discussed by Seiberg 
and Witten m- In many ways, it is similar to the 3- 
dimensional monopole plasma discussed by Polyakov [T2| . 

Unsal and Yaffe [T3] , using a double-trace deforma¬ 
tion of Yang-Mills at large N on x R^, argued that it 
prevents the spontaneous breaking of center symmetry. 
A similar trace deformation was used in the context of 
two-dimensional (confining) QED with unequal charges 
on X R |14j to analyze the nature of center symmetry 
and its spontaneous breaking. This construction was ex¬ 
tended to QCD with adjoint fermions by Unsal m, and 
by Unsal and others m to a class of deformed super- 
symmetric theories with soft supersymmetry breaking. 
While the setting includes a compactification on a small 
circle, with weak coupling and an exponentially small 
density of dyons, the minimum at the confining holon- 
omy value is induced by the repulsive interaction in the 
dyon-anti-dyon pairs (called bions by the authors). The 
supersymmetry is needed to calculate the contribution of 
the dyon-anti-dyon pairs , and, even more importantly, 
for the cancellation of the perturbative Gross-Pisarski- 
Yaffe-Weiss (GPYW) holonomy potential [17]. 

Shuryak and Sulejmanpasic [TS] have argued that in¬ 
duced by the “repulsive cores” in dyon-antidyon chan¬ 
nel also generates confinement, explaining it in a simple 
model. The first numerical study of the classical interac¬ 
tion of the dyons with anti-dyons has been recently car¬ 
ried in |22] . The streamline configurations were found by 
a gradient flow method, and their action assessed. This 
classical interaction will be included - for the first time 
- in our paper. 

Another major non-pertubarive phenomenon in QCD- 
like theories is spontaneous chiral symmetry breaking. 
Shuryak and Sulejmanpasic m have analyzed a num¬ 
ber of phenomena induced by the fermionic zero modes 
of the instanton-dyons such as the formation of clusters 
(molecules or bions) at high temperature and their col¬ 
lectivization, generating spontaneous breaking of chiral 
symmetry at low temperature. Faccioli and Shuryak [20j 
have started numerical simulations of the dyon ensem¬ 
ble with light fermions to understand the nature of the 
fermionic collectivization. We will provide an analytical 
analysis of these effects in the second paper of this series. 

Before we get into the details of the various approxi¬ 
mations to our analysis, let us try to provide some qual¬ 
itative answers to the three generic questions formulated 
above: (i) At T ^ T^, we will consider the action per 
dyon (anti-dyon) to be still large or S' ~ 4 whatever 
Nc, (ii) The dyon interactions will be of the order of 
ASint ~ 1 <C S. The quantum (one-loop) interactions 
are several times smaller and naively can be considered 
small. However they are quite non-trivial and the re¬ 
pulsion they provide would be our key finding, (iii) In 
general, the dyon plasma is strongly coupled and it is 
hard to treat it analytically. However we will argue be¬ 
low that in some window of temperatures (below Tc) one 


can still use the Debye-Huckel plasma theory. 

A major contribution to the understanding of the one- 
loop dyon interaction has been made by Diakonov and 
others la 0. They have found that at T > Tc their 
interaction with the surrounding QGP leads to a linear 
(confining) potential between the dyons, proportional to 
the perturbative Debye mass. Since in this work we will 
only consider the opposite case T < Tc, this will not be 
included in what follows. Key to the one-loop effect is 
the explicit quantum weight of the KvBLL instantons 
in terms of the collective coordinates of the constitutive 
dyons at all separations. The self-dual sector is charac¬ 
terized by a moduli space with a hyper-Kahler metric. Its 
volume element is given by the determinant of Coulomb¬ 
like matrix. We will refer to it as Diakonov determinant. 

In his first attempts to treat the dyonic plasma, Di¬ 
akonov kept only the one loop determinant, the volume 
of the moduli space, ignoring the QGP screening effects 
and - as we will discuss in detail - the even larger classi¬ 
cal dyon-anti-dyon interaction. Furthermore, he assumed 
that the attractive and repulsive terms induced by the 
determinant cancel out on average. We disagree on this 
conclusion as we detail below. Indeed, Bruckmann and 
others in m tried to generate configurations of randomly 
placed dyons using the determinantal measure, and ob¬ 
served that for the physically relevant dyonic densities, 
the determinantal measure develops negative eigenvalues. 
This makes no sense if the measure is to account for the 
volume of the dyonic moduli space. We will show that 
this issue may become resolved in a strongly correlated 
ensemble. 

It is well known that the separate treatment of self dual 
and anti-selfdual sectors is only justified in the context of 
supersymmetry where self-duality is dual to holomorphy. 
In QCD-like theories, the interaction between self dual 
and anti-selfdual sectors is strong and not factorizable. 
It is described semi-classically by a “streamline” with 
a classical inter-particle potential of order I/og, which 
is larger than the I-loop quantum induced potential of 
order a°. Furthermore, configurations with too strongly 
overlapping objects with small action, are not subject 
to the semiclassical treatment. To account for that one 
usually relies on the use of a “repulsive core” as in the 
instanton liquid model for instance. 

As we will discuss in detail, classical dyon-antidyon in¬ 
teraction |22j is about an order of magnitude stronger 
than the one-loop Coulomb effects. It generically leads 
to the dyon plasma in the strongly coupled regime, with 
^ 1 . yve will however focus on the very dense 
regime of such plasma, in which screening is strong 
enough that statistical mechanics of the ensemble can 
be treated by a variant of the Debye-Huckel mean field 
plasma theory. In such case the screening length is short 
enough to fence the system from strong coupling corre¬ 
lations and molecular-type instabilities induced by the 
streamline. The more dilute systems such as those ap¬ 
pearing at T > Tc, will not be discussed in this work, 
as they need more powerful many-body methods, such 
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as e.g. strongly coupled Coulomb plasmas many-body 
physics re-summations [231124j (and references therein). 
As we will show, in this case the free energy has a mini¬ 
mum at the “confining” holonomy value v = nT. 

In this paper we will detail the strongly coupled nature 
of the dyonic plasma. Our original results consist of (i) 
introducing the strong correlations between dyons and 
anti-dyons as described by the streamline [52]; (ii) show¬ 
ing that the determinantal interactions induced by the 
moduli space for dyons or anti-dyons are mostly repul¬ 
sive causing the moduli volume to vanish for randomly 
distributed dyons; (iii) showing that suitably organized 
dyons to account for screening correlations yield finite 
moduli volumes; (iv) deriving an explicit 3 dimensional 
effective action that account exactly for the screening of 
dyons and anti-dyons on the moduli space with strong 
inter-dyon-anti-dyon streamline interactions; (v) show¬ 
ing explicitly that the strongly coupled dyonic plasma 
is center symmetric and thus confining; (vi) deriving the 
Debye-Huckel corrections induced by the dyons and anti- 
dyons to the leading Pressure for the dyonic plasma and 
using it to asses the critical temperature for the SU(2) 
plasma; (vi) providing the explicit results for the gluon 
topological susceptibility and compressibility near the 
critical temperature in the center symmetric phase; (vii) 
deriving the scalar and charged structure factors of the 
dyonic plasma showing explicit screening of both elec¬ 
tric and magnetic charges at large distances with explicit 
predictions for the electric and magnetic masses; (viii) 
showing that the strongly coupled dyonic plasma sup¬ 
ports both electric and magnetic confinement. 

This paper is organized as follows: In section 2 we re¬ 
view the key elements of the dyon and anti-dyon measure 
derived in nig using the KvBLL instanton. The dyon- 
anti-dyon measure is then composed of the product of 
two measures with streamline interactions between the 
dyons and anti-dyons. We briefly detail the exact re¬ 
writing of the 3-dimensional grand-partition function in 
terms of a 3-dimensional effective theory in the SU(2) 
case. We also show that the ground state of this ef¬ 
fective theory is center symmetric. In sections 3-6 we 
show that in the linearized screening approximation the 
dyon-anti-dyon liquid still screens both electric and mag¬ 
netic charges, generates a linearly rising potential be¬ 
tween heavy charges and confines the large spatial Wilson 
loops. The t' Hooft loop in the dyon-anti-dyon ensemble 
is shown to be 1 modulo 0{as) self-energy corrections 
which are perimeter-like in section 7. Our conclusions 
are in section 8. 


II. INTERACTING DYON-ANTI-DYON 
ENSEMBLE 

A. The setting 

The first step is the introduction of the nonzero ex¬ 
pectation value of the 4-th component of the gauge field. 


which is gauge invariant since at finite temperature it en¬ 
ters the holonomy integral over the time period, known 
also as the Polyakov line. Working in a gauge in which 
(A4) belongs to the diagonal and traceless sub-algebra 
of Ac — 1 elements, one observes the standard Higgsing 
via the adjoint field. All gluons except the diagonal ones 
become massive. We will work with the simplest case of 
two color gauge theory Nc = 2, in which there is only one 
diagonal matrix and the VEV of the gauge field (holon¬ 
omy) is normalized as follows 

3 3 

(^ 4 ) = = 27rT:/y ( 1 ) 

where t^/ 2 is the only diagonal color generator of SU(2). 
At high T it is trivial with 1/ —>■ 0, and at low T < Tc it 
takes the confining value iz = 1/2. With this definition, 
the only dimensional quantity in the classical approxima¬ 
tion is the temperature T, while the quantum effects add 
to the running coupling and its A parameter. Since we 
are working near and below we will follow the lattice 
practice and we use the latter as our main unit. 

In the semi-classical approximation, the Yang-Mills 
partition function is assumed to be dominated by an in¬ 
teracting ensemble of dyons (anti-dyons) [H [5] . For large 
separations or a very dilute ensemble, the semi-classical 
interactions are mostly Coulombic, and are encoded in 
the collective or moduli space of the ensemble. For multi- 
dyons a plausible moduli space was argued starting from 
the KvBLL caloron |3] that has a number of pertinent 
symmetries, among which permutation symmetry, over¬ 
all charge neutrality, and clustering to KvBLL at high 
temperature. Since the underlying calorons are self-dual, 
the induced metric on the moduli space was shown to be 
hyper-Kahler. 

The SU(2) KvBLL instanton (anti-instanton) is com¬ 
posed of a pair of dyons labeled by L, M (anti-dyons 
by L,M) in the notations of |1|. Specifically M car¬ 
ries (-1-,-!-) and L carries (—,—) for (electric-magnetic) 
charges, with fractional topological charges Vm = v and 
vi = I — v respectively. Their corresponding actions are 
Sl = 2TTVmlas and Sm = ^TTVijas- 

The statistical measure for a correlated ensemble of 
dyons and anti-dyons is 
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( 2 ) 


The streamline interactions induced by the potential 
Vjjf) correlate the two otherwise statistically independent 
dyon and anti-dyon sectors. (Note that by the potential 
we mean the extra action and not the energy, thus no 
extra 1/T). Asymptotically, 
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(3) 


is a Coulomb-like classical interaction between dyons and 
anti-dyons. Here x^i and are the 3-dimensional coor¬ 
dinate of the i-dyon of m-kind and j-anti-dyon of n-kind. 
At shorter separations the streamline stops at certain dis¬ 
tance 0 ^ 1 />, we will refer to it as a “core size”. While the 
interaction is more complex than just electric Coulomb, 
it is proportional to the electric charges Q,Q. In general 
those are the (Cartan) roots of SU(Nc) supplemented by 
the affine root. They satisfy 

QmiQnj — ^m,nj-l ^m,n—l) (4) 

The dimensionality of G[x\ is {Ki + ... -\- and simi¬ 
larly for G[y\. Their explicit form can be found in 
In the SU(2) case there is only one electric charge. 

The semiclassical 3-density of all dyon species no = 

n-L + riM + ni+nM is 


no = 


dN 

d^x 


CT^, 


at 


(5) 


where C is a constant to be determined below (see ([5^). 
© can be re-written using the asymptotic freedom for¬ 
mula for SU(2) pure gauge theory with 27r/as(T) = 
(22/3) In (T/A) in terms of the scale parameter A. The 
dimensionless density 


We will use this fit as a basis for our running cou¬ 
pling. In particular, the action of the SU(2) caloron at Tc 
Sl+m{Tc) « 7.47 translates to the value of the coupling 
Cis(Tc) = 0.84. Since in this paper we only work in the 
confining regime of the holonomy with all dyon actions 
identical, the action per dyon is about 3.75. 

The repulsive linear interaction between unlike dyons 
(anti-dyons) found in acts as a linearly confining force 
in the center asymmetric phase, favoring the molecular 
or KvBLL configuration at T > Tc. This interaction 
stems from QGP thermal quanta scattering on the dyons. 
However, we will be interested in this paper in the center 
symmetric phase at T > T^, in which there is no QGP, 
we do not include this interaction. 

Since the classical ~ l/Os it dominates the quan¬ 
tum determinants, which include Coulomb interaction of 
order . On this point we differ from the argument pre¬ 
sented in [4] regarding the re-organization of (|^ in an 
extended quantum determinant. At large relative sep¬ 
arations between all particles the measure © is exact. 
It is also exact when each bunch of dyons or anti-dyons 
coalesce into a KvBLL instanton or anti-instanton at all 
separations. 

The above notwithstanding, the grand-partition func¬ 
tion associated with the measure ([^ 

— X f dy-Dpi^] (9) 

[K] 


np 

T3 ^ \t) 


( 6 ) 


is small at high T but increases as T decreases. With 
the exception of section [IH G[ where we will estimate the 
critical deconfinement temperature by including pertur¬ 
bative 0(a°)) effects in the dimensionless pressure, we 
will always assume the temperature to be small enough, 
so that the dyons effect are the dominant ones. The dyon 
fugacity / is 


/ 


np 

Stt 


(7) 


to order in the dyon density (see below). The 

absolute value of the parameter A appearing in the semi¬ 
classical formulae can be related to standard parameters 
like Ajj^, but this has no practical value since the ac¬ 
curacy with which they are known is too low to give an 
accurate value of the dyonic density. In practice it is 
obtained from the fit to the lattice instanton data per¬ 
formed in m in the range 0.5 < T/Tc < 3. The caloron 
action - the sum oi Sl and Sm ~ is then writen as 


describes a highly correlated ensemble of dyon-anti-dyons 
which is no longer integrable in the presence of the 
streamline. The case = 0 amounts to ^pp 

ZoZjj where each factor can be exactly re-written in 
terms of a 3-dimensional effective theory. We now ana¬ 
lyze ^ for the SU(2) case following and correcting the 
arguments in [1]. 


Kl Km Kj^ 

n n n 

[K] 2l,=l lM = l *L = 1 »M- = 1 

f fdPxLic fd^XMiM 

KlI Km'. Ki\ Km'. 
xdet(G[x])det(G[ 2 /]) e-^nnG-y) (iq) 


with G[a:] a {Kl + KmY matrix and G[y] a {Ki + Km)'^ 
matrix whose explicit form are given in mi]. 


B. Classical dyon-antidyon interactions 


Sl+m{T) = 


27r 


= —In f ^ ) 

3 [0.36 TcJ 


( 8 ) 


The explicit form of the Goulomb asymptotic in (10) 
for the SU(2) case is 
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\XM - VmI \^L - VlI \xm - Vzl \XL - VmI 

( 11 ) 


The strength of the Coulomb interaction in (111 is Cd /a^ 
and is of order l/ofs- It follows from the asymptotics of 
the streamline configuration. In Fig. we show the at¬ 
tractive potential for the SU(2) streamline configuration 
in the MM channel [52]. The solid curve is a numerical 
fit to the data given by 


= ( 12 ) 

with = — 1 in units of the critical temperature Tc 

and g set to 1 and ^4 = 30.9, B = 0.9072, C = 15.795. 
The dashed line corresponds to the Coulomb asymptotics 




MM 


(r) « - 


Cp 


OtsT 


(13) 


with Cp = A/Att = 2.46. We recall that in the un¬ 
combed DD potential, the asymptotic Coulomb inter¬ 
action corresponds to Cp = 2. The attraction in the 
streamline is stronger asymptotically owing to the rela¬ 
tive combing between the dyons. Fig. [T] shows that the 
DD core is about aj^p « 1/T. The second observa¬ 
tion is that one should not use the Coulomb asymptotic 
(the lower dashed curve) but the actual potential which 
correctly takes care of the dyons, as extended charged 
objects rather than point charges. 

Below the core value of app^ the streamline configu¬ 
ration annihilates into perturbative gluons making the 
parametrization (12) arbitrary. Throughout, we will 


parametrize the core by a constant, replacing (12) by 


^ppi’’’) = Spp {V{r)9{r — app) + V{app)9{app — r)) 

(14) 


this problem should be studied by methods more pow¬ 
erful than the usual mean field approximations, such as 
the Debye-Huckel theory. However, we will show below 
that when the dyonic densities are sufficiently large (and 
that implies the overall T of the ensemble to be suffi¬ 
ciently low), the screening mass gets large enough to put 
the effective - screened - interaction inside the domain 
in which the analytic Debye-Huckel theory becomes jus¬ 
tified. 

Furthermore, as we will detail below, the treatment of 
the repulsive core is in fact a rather sensitive issue. We 
chose the “most smooth” version of the potential, shown 
by the solid curve in Fig. Its Fourier transform pro¬ 
vides a smooth form factor in momentum space. We note 
that the actual streamline was only found for distances 
r > app « 1.2 (about A/v in the dyon units). The upper 
(blue) dashed curve is an example of an arbitrary param¬ 
eterization discussed in |22|, extending it to smaller val¬ 
ues of r. If one uses it, or even cut off the small r < app 
region completely - the approach known as hard core or 
excluded volume - the Fourier transform of the potential 
develops large oscillations. In this case the instability of 
the Debye-Huckel theory becomes stronger and its appli¬ 
cability domain sh rinks. 

The use of (12) in the repulsive channels ML and 
LAI approximates a smaller repulsion than Coulomb at 
shorter distances. A numerical investigation of these 
channels would be welcome. Note that both the mea¬ 
sure in (10) and the asymptotic (11) do not include the 
quantum corrections around the streamline configura¬ 
tion. Both of which should add more repulsion to the 
interaction between D and D. A leading quantum cor¬ 
rection to the asymptotic 0 follows by analogy from 
the Coulomb corrections emerging from the DD and DD 
determinantal interactions. In our case they are repulsive 
and amount to the shift 


Cp^Cp - -+0{al) (16) 

TT 

in the Coulomb constant. The relevance of this correction 
will be briefly discussed below. 


C. Qualitative effect of the one-loop moduli space 


’MM 


with 
-1-1 (repulsive). 
The ensemble 


= Si^l = -I (attractive) and = s^m = 


( 10 ) 


can be viewed as a 4-component 
dense and strongly coupled liquid. The quantity in the 
exponent, known as the classical plasma parameter 


Tdd — ^((^pp) 


Cp/agttpp 

3T^ 


(15) 


is not small. Its exponent can be even large. This 
implies that the “dyonic plasma” we want to study be¬ 
longs to a class of strongly coupled plasmas, with non- 
negligible correlations between the particles. So a priori. 


The volume element of the moduli space of self-dual 
SU(2) dyons is given by y/gnK = det C with gpK its 
associated Hyper-Kahler metric |2]. As we already men¬ 
tioned in the introduction, the one-loop determinant in 
the measure ^ must be positive definite for all config¬ 
urations. Furthermore, the positivity of all eigenvalues 
is required, since they have the meaning of the volume 
element in the corresponding subspace. As noted in m, 
this is not the case for ensembles with randomly placed 
dyons. These ensembles get denser and the positivity 
condition is only fulfilled for a very small fraction of the 
conhgurations. 
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FIG. 1: (Color online) Black solid line is the SU(2) DD (di¬ 
mensionless) potential versus the distance r (in units of 1/T). 
Upper (blue) dashed line is the parameterization proposed in 
Ref. [2^, the lower (red) (dashed) line is the Coulomb asymp¬ 
totics. 


In fact one of the main issues of the dyonic ensembles 
is the non-trivial character of the one-loop interaction 
induced by the Diakonov determinant. Before we show 
how this carries to our case through various fermioniza- 
tion and bosonization and diagrammatic re-summations, 
it is instructive to provide a qualitative understanding of 
the issues using simple explicit examples. 

Although it is well known, for completeness let us start 
with the simplest case of two dyons in the SU(2) theory 
with symmetric holonomy v = D = 1/2. Omitting the 
overall factors, Diakonov 2x2 matrix G reads 

/ 1 ± =F^— \ 

I VXi2 ' VXi2 \ 

G 2 X 2 N ^ (17) 

V 1±^ / 

\ ' VXi2 VXi2 / 

with xi 2 = 1 ^( 1 ) — ^( 2)1 the distance between the dyons 
in units of 1/v = I/ttT. The upper signs are for dif¬ 
ferent (ML) dyons, and the lower for similar (MM, LL) 
pairs. The metric-induced potential is thus V(xi 2 ) = 
—Indet G = —In (1 ± 2/(vxi2)) « =f2/(uxi2 ) is Coulomb¬ 
like at large distances. (At short distances the induced 
potential is proportional to ln(l/r) and not 1/r. There 
is no divergence in the partition function.) 

Let us now consider an ensemble of several (N = 8 ) 
dyons with JVm = A^l = 4 and set them randomly in a 
cube of size a. We then evaluate all inter-dyon distances 
and calculate detG[x] (which is now an 8 x 8 matrix) as 
a function of the Coulomb parameter e = l/( 7 raT). For 
each sampling, the determinant is a polynomial of e of de¬ 
gree JV. The results of 10 random samplings are displayed 
in Figj^by the dashed lines. For small e the determinant 
deviates from 1 in a non-uniform way. Some configura¬ 
tions are Coulomb attractive with detG > 1, while some 
others are repulsive with det G < 1 for small e. To first 
order, they average to zero for a large number of charges 
as there are equal number of positive and negative ones. 


At next order, the attraction is to win thanks to the gen¬ 
eral theorem of second order perturbation theory. How¬ 
ever, we observe that already for e = l/( 7 raT) ^ 0.2 the 
repulsive trend is dominant and det G < 0 for some sam¬ 
plings. This means that the moduli space of these con¬ 
figurations vanishes at the corresponding density. This 
sets an upper limit on the density of random ensembles 
of dyons 


n < uniax = 8 (0.2 7rTf ~ 1.98 (18) 

The lesson: Diakonov determinantal interaction for ran¬ 
domly placed dyons is strongly repulsive, reducing dra¬ 
matically the moduli space all the way to zero size for 
small e. It amounts to a strong effective core of order 

However this is not the end of the story. Let us look at 
the opposite case of a well ordered arrangement of dyons 
in the unit box. For that we pre-arrange the 8 dyons of 
the previous ensemble in a salt-like or fee configuration 
on the unit cube, and assess the corresponding det G. 
The result is shown in Fig. [^by the solid line. While the 
qualitative trend is the same - attraction at some interval 
of densities, changing to repulsion and then reaching zero 
at some density - the value of the maximal density to be 
reached is changed by a large factor of about 5^ = 125. 
Here is lesson number 2: the moduli space can be made 
much larger for the same inter-particle Coulomb strength 
e, if the correlations between charges are correctly taken 
into account. 

The overall lesson we get from those examples is the 
following: Diakonov’s original suggestion that attraction 
and repulsion would always cancel out is incorrect. Our 
analysis shows that ultimately the repulsion always wins 
and at some density the volume of the moduli space al¬ 
ways goes to zero. However, correctly implemented cor¬ 
relations between charges to maximize screening locally, 
can increase this critical density by about two orders of 
magnitude. 


D. Fermionization and Bosonization 


Following |4] each determinant in (10 1 can be fermion- 
ized using 4 pairs of ghost fields xl M’ Xl,m for the dyons 


and 4 pairs of ghost fields x^ XL,m foi' the anti-dyons. 
The ensuing Coulomb factors from the determinants are 
then bosonized using 4 boson fields vl,m,wl,m for the 
dyons and similarly for the anti-dyons. The result is a 
doubling of the 3-dimensional free actions obtained in [4] 


Sif[x,v,w\= J 

(|VxlP + |VxmP + Vul • Vwl + Vum • Vicm) + 


For the interaction part we note that the pair 

Coulomb interaction in (11) between the dyons and anti- 
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FIG. 2: (Color online) det G as a function of e = l/ijiaT). 
The dashed lines are for 8 dyons randomly placed in a cube of 
size a = 1. The solid line is for correlated dyons in a salt-like 
or fee configuration also in a unit cube. 


dyons can also be bosonized using the standard trick [T^ 
in terms of a and b fields. Here a and b are the un-Higgsed 
long range U(l) parts of the original magnetic field Fij 
and electric potential A 4 (modulo the holonomy) respec¬ 
tively. As a result each dyon species acquire additional 
fugacity factors such that 


L:e^+*'" M : Z : ( 20 ) 

These assignments are consistent with those suggested 
in HBUH] using different arguments. As a result there is 
an additional contribution to the free part (19) 


S2F[cr,b]=T J d^xd^y ( 21 ) 

X {b{x)V~'^{x - y)b{y) F a{x)V~^{x - y)a{y)) 
with V{r) defined in (12). The interaction part is now 


Si[v,w,b, a, x] 


d^x 


e + \xm - Xl\'^ +vm - vl) o’"" + 

(invi + Ixz, - XMp + vl - vm) + 

(47rv™ -h Ixm “ XhP + i ’ m - i ’ l ) -h 

^+b+^<yf 1 ^- _ ( 22 ) 

In terms of (|19l|22|) the dyon-anti-dyon partition function 
© can be exactly re-written as an interacting effective 
field theory in 3-dimensions, 


Z^jj[T]= J D[x]D[v]D[w]D[a]D[b]e-^^--^---^i23) 


In the absence of the screening fields a,b (23) reduces 
to the 3-dimensional effective field theory discussed in [4] 
which was found to be integrable. In the presence of ct, b 
the integrability is lost as the dyon-anti-dyon screening 
upsets the hyper-Kahler nature of the moduli space. We 
will investigate them by linearizing the screening effects 
in the symmetric state. 

Since the effective action in (23) is linear in the 
the latters are auxiliary fields that integrate 
into delta-function constraints. However and for conve¬ 
nience, it is best to shift away the b, a fields from ( 22 ) 
through 


wm — 6 + icr —t wm 

Wm -b-ia ^ Wm (24) 

which carries unit Jacobian and no anomalies, and re¬ 
cover them in the pertinent arguments of the delta func¬ 
tion constraints as 


-v-V^WM + 
47r 


= Tv ^(6 — ia) 

47r 


- — V^WL + = 0 

47r 


(25) 


and similarly for the anti-dyons. In [4] it was observed 
that the classical solutions to (25) can be used to inte¬ 


grate the w's in (23) to one loop. The resulting bosonic 


determinant was shown to cancel against the fermionic 
determinant after also integrating over the x^s in (231. 


This holds for our case as well. However, the presence of 
a, b makes the additional parts of (23) still very involved 
in 3 dimensions. 

After inserting the constraints in the 3-dimensional ef¬ 
fective action in (23), the ground state corresponds to 


constant fields because of translational invariance. Thus, 
the potential per unit 3-volume V 3 following from ([2^ af¬ 
ter the shifts (24) is 


- V/H 3 = 47r/ {vme 
+47r/ [vfne 


Wm-Wl 

WrOr-wr 


+ vie 


Vie ^ -rvi 

Wr-W^' 


(26) 


Note that if we did not perform the shift (24) then both 
the potential (26) and the constraints (25) depend on b 


and a making the extrema search for V more involved. 
Of course the results should be the same. For fixed 
holonomies Vm,i, the constant w's are real by ((25| as all 
right hand sides vanish, and the extrema of ( |26p occur 
for 


= ±\/V^ 

= Fyjvijvi, 




(27) 


(27) is only consistent with (25) if and only if vi = 
Vm = 1/2 and Vf = Vm = —1/2. That is for confining 



























holonomies or a center symmetric ground state. However 
and because of the constraint ( [ 2 ^ there is no effective 
potential for the holonomies in the interacting dyon-anti- 
dyon liquid. Indeed, by enforcing (251 before variation 
we have V/V 3 = —no, whatever the u's. On this point 
we differ from the arguments and corresponding results 
made in [ 3 ] where the constraints (25) were not enforced 
before the variational derivation of the holonomy poten¬ 
tial. Note that the alternative argument in [3] in favor of 
the holonomy potential fixes the number of dyon species 
Ki to be equal a priori, while (101 fixes it only on the 
average. 


III. LINEARIZED SCREENING 
APPROXIMATION IN CENTER SYMMETRIC 
STATE 

For the center symmetric ground state of the 3- 
dimensional effective theory, we may assess the correc¬ 
tion to the potential V to one-loop in the b, a fields. This 
is achieved by linearizing the constraints (251 around the 
ground state solutions. Specifically 


+ 2 /^ wm - ‘^fwL « 
-£ v 2 + 2 A WL-2fwM^0 


(28) 


and similarly for the anti-dyons. The one-loop correc¬ 
tion to V follows by inserting (28 1 in (23). The ensuing 


quadratic contributions before integrations are 


and the screening mass M = ^J^n^jT with | = 2 for 
SU(2). In Fig. we show the form factor (31) in dots 
line in units of T^. A simple parametrization is shown in 
solid line of the form 


g—p a£>£) 

V{p) « 4a-^— cos(pa£,s) 


(32) 


with a = TrColas and a core « l/T),. Inserting (32) 
into (30) and setting p = p/M yield 


InZiL/Pa = -V - 


M3 f (fp 

^ y (^ 


In 




with 


d{p) = ae “ddPQQ g (^MajjQp) 


(33) 


(34) 


The dominant contribution to the integral in (30) comes 
from the region p « 1 for which (34) can be approximated 


by d(l) = d. As a result (30) can be done approximately 


by fixing d and we have the classical contribution to the 
pressure 


^ = lnZiL/1^3 


m3 

nD + n{a) — 


(35) 


with 


= V - 47r/ 




(29) 


The coefficient of the b field appears tachyonic but is 
momentum dependent and vanishes at zero momentum. 


K,(d) = 


2 ^d 


1+f 


2-|d- 




1 -f 


-4 


(36) 


(36) is seen to vanish for d = 0 or in the absence of DD 
interactions. Near T^, the screening mass is M « ueIT'c 
(see below), thus 


A. Pressure 


Carrying the Gaussian integration in 6 , cr in (29) yields 
to one-loop 


InZiL/Hs = -V - - 


d^p 

(^ 


In 


1 - 


V^{p) p^M^ 


16 (p2 -H M^y 


with V (p) the Fourier transform of (12) 


d = {nCD/as) cos {MaoD) « -0-52 (37) 

For |d| < 4 the 1-loop contribution to the pressure from 
the charged DD dyons is real with no dimer or molecular 
instability. The large core produced by the form factor 
(34) is considerably screened by the large dyon density as 


(30) 


captured by the large dielectric constant 1 /k(— 0.52) « 
5.26 in 

The correction in (35) to the free contribution is a 
Debye-Huckel correction [23] (and references therein). A 
simple but physical way to understand it is to note that a 
screened Coulomb charge carries a lower constant energy 




dr sinrVj^^{r/p) 


(31) 




1 


Ait\x\ At:\x\ 


M 

Att 


(38) 
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FIG. 3: (Color online) The dots show the fo rm factor, the 
ratio V(p) ■ (p^/47r) of the Fourier transform of (121 to that of 
a pure Coulomb law versus p/T. The thin line is its parame¬ 
terization. See text. 


The Debye-Huckel as a mean-field estimate for the pres¬ 
sure follows 



. ' ■ . 

0.6 0.7 0.8 0.9 1.0 1.1 


Vdh _ udM _ M3 

T 47rr Stt 127r 


(39) 


where no = M^r/2 is the density of charged particles 
(see below). The standard Debye-Huckel limiting result 
for a multi-component ionic plasma in 3 spatial dimen¬ 
sions is shown on the right-most side of (39). 

The correction in (35) is considerably reduced by the 
large screening through the effective dielectric constant 
played by 1/^(5) « 32/(15d^) for d ^ 1. In particular 
!/«;(—0.52) « 5.26 3> 1 as noted earlier. It can be recast 
in the form 


^ci ~ K(d) 4 
—r = nn H-;= np. 


(40) 


with fiD = ud/T^- Using Ma^p, « cte/T^ ~ 1/(0.71)^ 
for SU(2) we have fip « 1, so that Vc\/T'^ « (1 -I- 0.01)). 
The screening corrections are small of the order of 1% 
thanks to the large dyonic densities. 

The limitations of the Debye-Huckel approximation are 
readily seen from (30). In Fig. |4k we plot the argument 


of the logarithm in the last term of (30). The differ¬ 


ent curves from top to bottom follow from Mapp = 
1.5,1,0.7, 0.56 respectively. The smaller the Debye mass 
M the stronger the dip. For Mapp < 0.56, the argu¬ 
ment of the logarithm becomes negative resulting into an 
ITT contribution to the pressure and thus an instability. 
This is a clear indication of a well known phenomenon: 
the Debye-Huckel approximation is in general inapplica¬ 
ble for strongly coupled plasmas, and the interaction me¬ 
diated by the streamline is strong. Only a large enough 
density of dyons, producing sufficiently strong screening, 
allows for the use of the Debye-Huckel theory. In Fig. 


FIG. 4: (Color online) (a) The argument of the logarithm 
in the last term of (301 versus the dimensionless momen¬ 
tum p, for different values of the dimensionless Debye mass 
Mapp) = 1.5,1,0.7,0.56, top to bottom. As the screening 
mass decreases to its critical value, the lower (green) curve 
touches zero. The smaller values of M leads to a negative 
argument of the logarithm, thus an instability, (b) A semi- 
logarithmic plot of the integral entering in ( |30[ ) as defined in 
(41) as a function of Mapp,. The decrease is steady from 
its maximum at the critical value of the screening mass or 
Mappi = 0.56. 


we show how the total integrated contribution to the free 
energy changes as a function of the dimensionless Debye 
mass Mqjpp) 


{Mappp,)^ 



V^{pM) 

16 (p3 -I- 1)4 


(41) 


The main lesson is that beyond the critical value of the 
screening, this contribution becomes rapidly very small. 
This is consistent with the analytical estimate above. 
This justifies the use of the Debye-Huckel mean-field 
analysis in general, and the use of the semi-classical ex¬ 
pansion in particular. 


B. Beyond the Debye-Huckel theory 

The unravelling of the Debye-Huckel approximation 
may be due to corrections to an interacting Coulomb sys¬ 
tem, such as 1/ core corrections; 2/ dimer, tetramer and 
so on many-body interactions. The large core corrections 
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were already identified and discussed above and yield a 
substantial reduction in the Debye-Huckel contribution 
near the critical value oi Ma 


DD 


0.56. 


Bound state corrections in the form of electrically 
charged LL or MM dimers, or electrically neutral 
LLMM tetramers commonly referred to as instanton- 
anti-instanton m olecule s, can bind through the stream¬ 
line interaction ( TT][l2 ). The combinations LM and 
ML are repulsive. The binding energy in a dimer is 
^DD ~ {Cd/'^)/{ oisaDo) = T. The dimer enhancement 
is expected to be of order ^ ^(Co/3)/as ~ 2.72 for 

T f^Tc using the reduced effective Coulomb coupling. As 
we noted earlier, this enhancement becomes substantially 
larger at high temperature as as decreases with the onset 
of dimerization set at about Og^crit = 1^)14: « 0.67. 

At this coupling which occurs above T^, the Coulomb 
dimer enhancement factor is « 3.57. 

In sum, the dyons and anti-dyons form a Coulomb liq¬ 
uid with strong short range correlations induced by both 
the finite cores and bindings. The liquid supports center 
symmetry and confines. The deconfinement transition 
is characterized by clustering into charged dimers and 
possibly uncharged and topologically neutral tetramers, 
forming mixtures with the restoration of center symme¬ 
try. Coulomb mixtures present rich phase diagrams |26j . 


C. Dyonic densities 


(30) can be readily used to assess the dyon densities 


Km and (and similarly for Km and Ki) in the center 
symmetric vacuum with screening dyons-anti-dyons. For 
that we need to change / —)■ '/JmJl and take derivatives 
of (301 with respect to In/^.L separately and then setting 
them equal by bulk charge neutrality. The result per 
species is 


... 1 

/r = 3»D + »(a)5^ 


(42) 


for all dyon and anti-dyon species. 

Each dyon (anti-dyon) is characterized by an SU(2) 
core of size p « l/(27rT v) « 0.33 fm in the center sym¬ 
metric phase with v = \/2 sX T = 1/fm. The De¬ 
bye length Xd = 1/M « ^T/2nD ~ 0.70 fm is about 
twice the core size. The classical Coulomb ratio for the 
DD, DD pairs with a core of 2p is about 


DD.DD = 


2tt{2p)T 2 


1 
' 4 

1 . 


(43) 


is small. Recall from (15) that T^,^ « 1. The Coulomb 
DD, DD interactions are quantum and of order with 
strength I/tt as can be seen by expanding the exponen¬ 
tial form of the determinantal interaction in (17). The 


dyon-anti-dyon ensemble is close to a strongly coupled 
4-component Coulomb liquid. Since the measure for the 



FIG. 5: Topological susceptibility in units of T versus T/T^ 


unlike dyons is exact, it is valid even in the dense con¬ 
figuration. It is only asymptotically exact for like dyons. 
For the dyons and anti-dyons the streamline is numeri¬ 
cally exact at all separations outside its core. However its 
corresponding quantum determinant was not calculated. 
Only a qualitative correction was argued in (24). 


D. Gluon condensates and susceptibilities 

The topological charge fluctuates locally in this dyon- 
anti-dyon model. The topological susceptibility at 1-loop 
follows from (30) through the substitution / —>• /cos(0/2) 
both in V and also M ^ M •\/cos(0/2). At finite vacuum 
angle 0 and in leading order we have 


FF) =- 


T dlnZiL 


0 E, 89 


/1 M^T _\ 

sin(6»/2) ) 


(44) 


Thus the topologicl susceptibility 




-) {tidT) (45) 


in leading order. Since the dyons carry half the topolog¬ 
ical charge (45) shows that the topological fluctuations 
are Poissonian to order The behavior of Xt/T'^ 

versus T/Tc is shown in Fig. with no defined in (56) 
below. 

The gluon condensate to 1-loop in the screening ap¬ 
proximation follows from 


1 


(^^)o - - 


167r2 ' 4 0 

_ _T_ f2^ 

Fk V a. 


T dlnZiL 
47rVi3 81/as 

K{a)M^ \ 
Stt J 


no 


(46) 
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FIG. 6: The electric (solid-black), magnetic 

(dashed-blue) and (Euclidean) energy density 
(dot-dashed-brown) in units of T versus T/T^. See text. 


which is non-Poissonian because of the scale anomaly. 
The compressibility of the ground state is 


E. Electric and magnetic screening masses 

The center symmetric phase of the dyon-anti-dyon liq¬ 
uid screens the long-range U(l) gauge fields left un- 
Higgsed by the holonomy Ai{oo)/2'KT = vT^/2. The 
electric and magnetic correlations in these Abelian U(l) 
charges can be obtained by introducing the U(l) Abelian 
sources rim,e for the magnetic and electric charge densi¬ 
ties 


z{x) — 'y ' Qn 


- Xi) 


(50) 


with \Qm,e\ = 1, and shifting the U(l) fields cr —cr -I- 
and 6 —>■ 6 -b in the 3-dimensionsl effective action. To 
1 -loop the generating functional for the charge density 
correlators is 


CTv = 


V 

T \ Vl67r2y 


ZiL[rimM= J 


(47) 


T 

IGtt^ 


no 


kM^\ 
Stt / 


no 


3kM3\ 
IOtt ) 


for the connected correlator. 


We can use (35) and (46) to extract the electric 
and magnetic (^)g condensates in the dyonic ensemble. 
For that we note that the energy per volume in Euclidean 
space follows from (35) through 


which is Gaussian in the sources and therefore readily 
integrated out. Thus 


\nZiL[iim,Ve] = - j 7 ^^ (.52) 

V / o*-c 


with the electric and magnetic density correlators follow¬ 
ing by variation, 




The results are 


47rr 

47rT 


2 a' /I TT 

+3-(l + -#) 


ai 




at 


no 


no 


(48) 


kM^\ 
Stt ) 
kM^\ 
Stt ) 
(49) 


with a'g = das/d\r{T. In Fig. we show the behav¬ 
ior of the chromo-electric condensate (solid-black), 

the chromo-magnetic condensate (dashed-blue) and 
the (Euclidean) energy density (dot-dashed- 

brown) in units of T versus T/T^ in the center symmet¬ 
ric phase. We used the dyon density fixed in ( [^ below. 
The chromo-magnetic condensate is about constant in 
the range of 0.6 < T/Tc < 1 while the chromo-electric 
condensate decreases monotoneously. The condensates 
are about equal and opposite near a point supported 
by the lattice extracted condensates in . 


1 / 

ue(p) = ^ {\PmAp)\‘ 


TMA 


2^4 


4(p2 + M2)2±aMV 

(53) 


The upper sign is for magnetic and the lower sign for elec¬ 
tric. In x-space, (53) can be inverted by Fourier trans¬ 
forms. The result for the electric correlator in spatial 
coordinates is 


TM^ 


cos 



-bsin I M\x 


( 1 -f) 


(54) 


The magnetic correlator follows by analytical continu¬ 
ation through the substitution d —?► —a in (54). The 


electric screening masses Mm,e follow from the large dis¬ 
tance asymptotics. Using our estimate of d « —0.52 < 0 
from the Debye-Huckel analysis above, we have 
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Me 

~W 



Mm 

~W 



« 0.93 


(55) 


with = 2 nD/T. Also the arguments below show that 
M = ue/T. Combining these two results allow us to fix 
C in ® above. Indeed, at T^. the SU(2) lattice results 
give ^foE ~ 0.71. So ^ now reads 


^ o O^sjTc) 
T3 ^ a^iT) 


3(Tc) 


s{T) 


(56) 


which gives Me ~ 1.4 Tc and Mm « 1.8 Tc, both of which 
are remarkably close to the reported SU(2) lattice results 
in the vicinity of the critical temperature m- In Fig. [7] 
we display the results (551 for Me,m/T in the range 
(0.5 — 1) Tc versus T/T^. The points at T > Tc are shown 
for comparison. We note that the electric mass drops 
down at Tc. In the region we study Mm > Me, while 
above Tc, in a more familiar QGP region, Mm < Me- 
This switching of the magnitude of the two screening 
masses is better documented in lattice works with the 
SU{3) gauge group. It has a simple explanation in our 
case. Since at T > Tc the dyon density drops it follows 
that M decreases as well. As a result, the form factor in 
Fig. [^is probed at smaller momentum p « M (larger dis¬ 
tances) making a{p ~ M) in (34) switch from negative 
(T < Tc) to positive (T > Tc). From (53) it follows that 
the expressions for Me,m in ( |55[ ) are now switched with 
Mm lighter than Me- A simple estimate of the critical 
temperature at the crossing follows from the vanishing 
of (37) or MacD = '^/2- This translates to a critical 
dyon density « n'^T^/S which is consistent with our 
estimate of Tc below (see Q)- 

Finally, we note that the value of as{Tc) « 0.84 ex¬ 
tracted from the cooled caloron data in m is also consis¬ 
tent with the reported value from bulk thermodynamics 
in |28j . In the dyon-anti-dyon Coulomb liquid the corre¬ 
lators are modified at intermediate distances as we now 
detail in terms of the static structure factors. 



FIG. 7: (Color online) The electric Me/T (dashed line) and 
magnetic Mm/T (solid line) screening masses in (55) versus 
T/Tc- The points are SU(2) lattice data from l27l shown 
for comparison, (blue) circles are Me/T, (red) squares are 
Mm/T- 



GE(p) 

GM(p) 


FIG. 8: (Color online) The electric and magnetic structure 
factors (|53|) as a function of p/M- 


Gm,e{p) 


Gm,e{p) ^ P'^ _ 

riE /2 (p^ -I- 1)^ ± dp^ 


(58) 


with p = p/M- We note that the pre-factor in ( [^ in¬ 
volves two static electric or magnetic exchanges with an 
identical screening mass M. The charged structure fac¬ 
tors vanish as Gm,e{p) ~ P^- For large momenta or 
p ^ 1 both structure factors asymptote one from below 
as shown in Fig. The magnetic hole is slightly smaller 
than the electric one around the same pairs. The absence 
of oscillations in the structure factor, is a consequence of 
our linearized approximation. 

To characterize further the 4-component plasma of 
dyons and anti-dyons we define the scalar static pair cor¬ 
relation function 


F. Static structure factors 


The charged structure factor between pair of magnetic 
or electric charges is (53) which can be re-written as 


Gm,e{p) 



Nrr,,. 

E 

i=i 


Qr\ 


ik-: 


(57) 


Thus 


Gs(a;) 



N 




Xi - 



(59) 


normalized to the total number of particles N. (59) de¬ 
fines the probability to find two particles a distance |a;| 
apart. Its Fourier transform 


Gs{p) 


1 

N 


N 


i=i 



(60) 
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is the scalar structure factor. 


(60) can be evaluated by switching / —>• f + df{xi) 
in (10) and then linearizing the resulting effective action 


around the mean-density. Specifically, we can re-write 
the linearized constraint (281 formally as 


{WM - Wl) = 


1 




Ag -I- 4(5/ \ 47r 


{b — ia) (61) 


with Ag = —TV^/47r -|- 4/ and use perturbation theory 
to expand the denominator in (61) to order 0{6f^). The 


result can be formally written as 


/ rl^k 

j^G{p,k){b - ia){k) 


(62) 


Inserting (621 into the potential (261 yields a quadratic 


action in b and a. Integrating over the latters yields the 
1-loop determinant or the effective action for S = Sf/f. 
Specifically 


det(l + S[5]) = e™"(i+«[^l)«eT^®[^l 


(63) 


with the quadratic effective action for the scalar fluctua¬ 
tions as 


TrS[5] = - 


with to order 0 {a^) 


d^p 

{2tt) 


3 (5(p)G ^{p)S{-p) 


(64) 


G-^P) 


ud +AciM° 


d^k 

(^ 

1 


(fc2-pM2)4((fc-Pp)2-HM2)2 
2fc2 1 


(P +M2)3 ((fc-pp)2 +M2)2 


(65) 


The njj contribution in (65) follows from the expansion 


of the leading contribution V using arguments similar 
to those used for the derivation of the dyonic densities 
above. 


The scalar structure factor follows from (65) through 
the normalization 


Gsip) = ^ (l<5(p)|^) = no G{p) 


( 66 ) 


We note that the small momentum fluctuations in 6 f 
couple to the sound-like modes. Specifically, 


is dominated by a massless pole at zero momentum with 


d « l-k8d^ 


(M\ 

r d^k 

r d 2 k^ 1 


1 (27r)3 

_(A:2 + 1)6 (fc2 + l)5_ 


Alternatively from the pressure (35) we expect 




dVci 

Tdn 


1-f 


D 


K{a) 

47r 


( 68 ) 


(69) 


with k(( 5) « 15 5^/32 in leading order and in total agree¬ 
ment with (68). Also, at large momentum (66) asymp¬ 
totes 65(00 )A 1. The slight super-luminal character of 
( [6^ reflects on the fact that dyons are in essence Eu¬ 
clidean configurations with no physical particle realiza¬ 
tion. 


G. Estimate of the critical Tc 

The total thermodynamical pressure of the dyon-anti- 
dyon liquid consists of the classical and non-perturbative 
contribution (40) plus the perturbative holonomy poten¬ 


tial known as Gross-Pisarski-Yaffe-Weiss (GPYW) po¬ 
tential [m, plus the purely perturbative black-body con¬ 
tribution (ignoring the higher order 0 {as) quantum cor¬ 
rections). Specifically (A^ = 2) 


'Ptot 

2^4 


' no 

TT^ 


K{a) 


n 


SttV2 

1 


D 


Ni- 




(70) 


The Debye-Huckel contribution is of order A^, while the 
leading classical contribution is of order N'^. So screen¬ 
ing and large Ac are not commutative. For the SU(2) 
case of interest, the transition temperature from the 
disordered phase {v = 1/2) to the ordered phase {v = 0) 
occurs when the first three contributions in (70) cancel 
out. Thus 


no 


K^Ot) 

■i-Ky/2 




12 


(71) 


For k(-0.52) 


^tHyi 


0 . 88 . 


0.19, the critical density is 
Since the SU(2) electric string ten¬ 


sion is aE = TM = \j2nDT (see below), it follows that 
= 6i/v/r « 0.88 which is somehow larger than 
the SU(2) lattice result T^j= 0.71 |29j. 


IV. POLYAKOV LINES 


Gs(p) 


C2p2 


(67) 


To probe the confining nature of the dyon-anti-dyon 
liquid in the 3-dimensional effective theory we will 
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compute explicitly the expectation of a heavy quark 
through the traced Polyakov line and the corelator of 
a heavy quark-anti-quark pair through the correlator of 
the traced Polyakov line and its conjugate at fixed spa¬ 
tial separation. The insertion of these charges in the 
dyon-anti-dyon liquid modifies the ground state through 
solitonic solutions around these sources. 

In this section we present a new derivation of the per¬ 
tinent solitonic equations for the SU(2) case that makes 
explicit use of the presence of the long range U(l) 6 and 
a helds. In the linearized screening approximation, we 
show that the solitonic equations for the heavy source 
probes are in agreement with those established in ^ us¬ 
ing different arguments. 


setting b = a = 0. (75) is a Poisson-Boltzmann type 
equation. It is also referred to as an elliptic and periodic 
Toda lattice mm- The solution is a local Debye-like 
cloud around the inserted heavy quark 


This causes almost no change in the vacuum holonomies 
Vm,i- Thus, after the shift 


(L(xi)) « =0 


(77) 


A. (L) 

In the SU(2) case the Polyakov line consists of inserting 
a heavy quark whose free energy consists of its Coulomb 
interactions with all the Coulomb charged dyons and 
anti-dyons. Specifically, the traced Polyakov line before 
averaging is 


L(xi) = 


(72) 


27r ifJ,M+ 2T X]i(^ 1^1 I ) 


with fiL — fJ-M = Vm- When averaging using the ensemble 
@ it is clear that each of the contributors to the string of 
factors in (72) will match its analogue from the measure 


and re-exponentiate. For instance the first contribution 


in (72) re-exponentiates through the substitution 




2T\xi-x 


(73) 


The extra Coulomb factors can be re-defined away by 
shifting 


b-^b + 


2 r|a;i — x\ 


(74) 


thereby changing the constraint equation (25) to 


Att 


^) = 0 (75) 


and similarly for the second contribution in (721 with 


L o M. The effect of the first contribution in the 


Polyakov line (721 is to add a source term to the con¬ 


straint equation for wm- It is in agreement with [1] after 


B. (LL^) 

The preceding analysis can also be applied to the cor¬ 
relator of two heavy quarks through LL^ which consists 
now of 4 contributions before averaging 

L{xi)L\x2) = 

m,n—M,L 

(78) 

with the pertinent Coulomb free energies Fm( 2 ^ 1 , 2 ) fol¬ 
lowing from ( |72[ ). When averaged over the measure (§, 
each of the factors in (78) can be matched with its ana¬ 


logue in the measure. The preceding observations show 
that the Coulomb factors in the probing correlator can 
be paired with 




2T\x-xi\ 2T|x-X2 I 


(79) 


A rerun of the preceding arguments shows that the con¬ 
straint equations acquire now two source contributions, 
one for each of the heavy quark inserted 

47r 

= ^V^(6 - ia) ^[5^{x - cci) - 5'^{x - X2)\ 

47r z 

-^V2u;L + /(e“'^-"'"-6“'“-“'^) =0 (80) 

47r 

Since 'S/‘^\/\x — X 2 \ = —4TrS^{x — X 2 ), we can symmetrize 
(80) by shifting S^{x — X 2 ) from the first to the second 


equation through 


Wm,L 


Wm,L 


t/T 


2 |a; - X2\ 


(81) 


with unit Jacobian. The symmetrized ( 80||81[ ) equations 
are in agreement with those established in [4] for the 
SU(2) case after setting b = a = 0. In this case the 
solution is peaked around the heavy quark sources 



























15 


{WM - Wl){x) 


27ri 

~Y 


d^p 

(27r)3 p2 _j_ ]\^2 


<82) 


Inserting this back in the expectation value of the corre¬ 
lator (78) yields asymptotically 


{L{xi)L\x 2 )) « 


(83) 


in the 3-dimensional effective theory in agreement with 
the result in [3] . In 4-dimensions (831 translates to con¬ 
finement of the electric charges with the electric string 
tension cr^ = MT. The additional Coulomb screening 
in (§ does not affect the asymptotics of the linearly ris¬ 
ing heavy quark potential to leading order. The dyon- 
anti-dyon Coulomb liquid still electrically confines in the 
center symmetric phase. 


long range U(l) a and b fields in the leading order ap¬ 
proximation. Thus (TrkC(C)) « q-<^m S jg saturated 
by the pinned soliton, with the magnetic string tension 
om = (Je = MT. This result is expected from the equal¬ 
ity of the electric and magnetic masses in (55). 

A simple understanding of this result is as follows: 
while a heavy quark sources an electric field, a large spa¬ 
tial Wilson loop sources a magnetic field by Ampere's law 
which is classically composed of all the magnetic poles 
fluxing S as is explicit in (85). The typical contribution 
to (84) for a planar surface in the xy-plane is then 


dIsB-dS 


Is{BAx,v)B,i0,0))dS 


( 88 ) 


by keeping only the first cumulant in the average and 
using translational invariance for large S. In this limit, 
S acts as a uniformly charged magnetic sheet made of 
magnetic dyons classically, so that 


V. SINGLE-WINDING WILSON LOOP 

To study the large spatial Wilson loops we use the same 
observations made above in the presence of the U(l) fields 
a and b. As an observable the traced spatial Wilson loop 
of area S supported by the spatial contour dS = C reads 


TrW(C') = e^dsB+-dS ^iS^B. dS 
and sources the static magnetic field 


(84) 


- ±2^ Vz I |3 


(85) 


When averaged using (§, the spatial Wilson loop ( [M| 
modifies the additional U(l) fugacity factors in the dyon 
sector. Their contribution follows again by shifting b =F 
itr —>■ 6 =F i((T — 77±) in the constraint equations with 


'q±{x) = ± / dSy 


x-y 

2 \x-y\ 


( 86 ) 


As a result, (25) in the presence of (84) are now modified 
to read 


T . 




= —V^p+ix) 


An 


47r 


0 = 0 


(87) 


for the first contribution and similarly for the second 
contribution in (84) with 77+ — p-. After choosing 


the spatial Wilson loop to lay in the x-y plane through 
V^77± = ±An5’{z), the results (87) are in agreement with 
those derived in a for the SU(^ case but without the 


{B,{x, y)BM 0)) « ^ ^ ^ (89) 

where the variance in the magnetic charge is assumed 
Poissonian. The magnetic charge density per unit 4- 
volume is (TM)^/2. The typical magnetic charge per 
unit area is then about its square root or {Qm) /S « TM. 
Thus 


(e^SsB-dS^^^-kMTS (go) 

which is the expected behavior up to a factor of order 
one in the string tension. 


VI. DOUBLE-WINDING WILSON LOOP 


Recently it was pointed out in m that a co-planar 
and double winding Wilson loop in the SU(2) pure gauge 
theory version of the model discussed by Diakonov and 
Petrov [3] shows an exponentl fall-off with the sum of 
the areas. In contrast lattice SU(2) simulations appear 
to show an exponential fall off with the difference of the 
areas. The main observation in m is that the solitonic 
configuration contributing to the single-winding spatial 
Wilson loop as for instance from our linearized version 
with & = CT = 0 in (87), factors out in the the double¬ 


winding and co-planar Wilson loop. 

For two identical loops with Ci = C 2 = C, we have 
the formal SU(2) identities [32] (and references therein) 


(Tr W(C)f = Trs {W{C)) + Tta (VF(C')) 
Tr(VF(C)2)= Tts{W{C))-Ma{W{C)) (91) 


The simple trace Tr is carried over the fundamental rep¬ 
resentation of N-ality fc = 1 as in (84), and Trg^^ are 
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carried over the symmetric and anti-symmetric of N-ality 
k = 2 (modulo 2) representations of SU(2) respectively. 
The identities (921 are commensurate with the Young- 
Tableau decomposition. In the dyonic plasma considered 
here, the k-string tensions ak in the linearized plasma 
approximation are identical to those derived in [4] with 
^kl^i = sinA:7r/2 for SU(2) with ai = aE- For k = 2 we 
have cr 2 = 0 and the second identity in (92) implies for 
large loops 


(Tr(lT(C)2)) = (Trs(IF(C)))-l (92) 

We have set all self-energies to zero for simplicity as they 
depend on the subtraction procedure. ( [^ is consistent 
with the doubly traced Wilson loop as dominated by the 
k = 2 modulo 2 colorless di-quark-like {qq) or baryon-like 
configuration in SU(2). In the dyonic plasma, the dou¬ 
ble Wilson loop with Ci = (72 is a bound colorless state 
with zero-size that is strongly correlated within the dyons 
cores and therefore is consistent with the arguments pre¬ 
sented in m- 

For largely separated loops Ci ^2 of arbitrary sizes but 
still lying in the spatial directions, clearly 


{Tt(W{Ci)W[C2 ))) « (93) 

for {Ai + A 2 ) < ^12 where Ai ,2 are the planar ar¬ 
eas supported by Ci ^2 separately, and A 12 is the min¬ 
imal area with boundaries Ci and C 2 - The main is¬ 
sue is what happens for the same doubly wound SU(2) 
spatial Wilson loops when A 12 < {Ai A 2 ) 1 Here 
we note that LL and MM dimers carrying (—2,0) 
and (+2,0) (electric, magnetic) charge assignments could 
cluster around the probe qq (baryon) and qq (anti- 
baryon) configurations respectively, to form neutral 
molecular bound states with masses that scale with A 12 
instead of (Hi + H 2 ). They are commensurate with the 
massive off-diagonal and charged gluons Higgsed by the 
holonomy and dropped in the dyon liquid analysis. These 
configurations were not retained in [3]. 


flux across the loop. The t' Hooft loop corresponds to 
a dual charge in the center of the gauge group running 
around a closed loop and measuring the chromo-electric 
flux across the loop. The i' Hooft loop is the dual of the 
Wilson loop. 

In the temperature range 0.5 Tc < T <Tc confinement 
is still at work and we expect the temporal Wilson and t' 
Hooft loops to exhibit behaviors similar to those in the 
vacuum state. In section |V] we have explicitly checked 
that the closed spatial Wilson loop obeys an area law. 
The temporal Wilson loop is not amenable to our dimen¬ 
sionally reduced and Euclideanized effective field theory. 

The t' Hooft loop V(C) enforces a gauge transforma¬ 
tion ric which is singular on a closed curve C. If a curve 
C' winds ncc times around C then 


{C)W{C')V{C) = e*^"<"c'IF((7') (94) 

V{C) amounts to a multi-valued gauge transformation 
on the loop (7, 


Hc(6» = 27r) = (95) 

with 6 an affine parameter along C. A simple choice is 


Hc(x) = 


(96) 


where <pc(cc) is a multi-valued scalar potential for the 
magnetic field Be generated by a loop of current jc run¬ 
ning along (7, and Q = (1,1,..., —N^ + 1) a Cartan gen¬ 
erator of SU(Nc). An alternative construction using a 
discontinuous solid angle was discussed in [37l [38] . The 
effects of (96) on an Abelianized Wilson loop is 


ds-A'^ ^ ds-{A-I^QBc) ( 97 ) 

with Bq = —Vipc- Note that since ipc is multivalued we 
have V x Be = 47rjq. If we normalize the loop current 
je such that 


VII. T' HOOFT LOOPS 

In an important study of the nature of confinement 
in gauge theories, t' Hooft [S3] has introduced the con¬ 
cept of a disorder operator or t' Hooft loop to quantify 
confinement in the Hilbert space of gauge configurations. 
The t' Hooft loop is a canonical operator much like the 
Wilson loop. In a Lorentz invariant confining vacuum, t' 
Hooft has argued that the temporal Wilson loop and the 
t' Hooft loop cannot exhibit an area law simultaneously. 
The temporal Wilson loop obeys an area law, while the 
t' Hooft loop obeys a perimeter law. 

Physically, the Wilson loop corresponds to a color 
charge in the fundamental representation running around 
a closed loop and measuring the the chromo-magnetic 


ds ■ Be = 47r 


ic 


lA{C') 


dS ■ je = -nee 


(98) 


then (97) reduces to 


In the space of gauge configurations, the gauge transfor¬ 
mation Qe is inforced through 


V{C) = d^xTT{EiDi{Q<fic)) (200) 

For SU(2) we have 
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V{C) = ^ (^01) 


where the latter substitution E —>■ iE follows in Eu¬ 
clidean space. With this in mind, the expectation value 
of the t' Hooft loop in the dyonic ensemble involves a 
string of sources to be inserted in (101. In leading order 


N+N 


V{C)-. n 


2 = 1 
N-\-N 


( 102 ) 


i=l 

Thus {V{C)) = 1 modulo 0(as) Coulomb-like self-energy 
corrections which are perimeter-like in general. 

Finally, the Polyakov line as a Wilson loop around 
the periodic temporal direction has a dual Polyakov loop 
with a dual magnetic charge in the center. In the con¬ 
fined phase, the temporal component of the gauge field 
asymptotes fixed electric-type holonomies, while its 
dual A 4 asymptotes zero dual magnetic-type holomies 


thanks to parity. A rerun of the arguments in IV A shows 
that while (T(x)) = 0 in (771 as expected in an Euclidean 


and confining thermal state, its dual does not vanish, i.e. 


(|l(x) = Tr ■I'o = 1 (103) 

again modulo 0{as) Coulomb corrections. This behavior 
is consistent with the one reported on the lattice for Nc = 
2,3 [21]. 


VIII. CONCLUSIONS 

The central theme in this paper is non-perturbative 
gauge theory for temperatures in the range (0.5 — 1) Tc 
modeled by a dense plasma of instanton-dyons. The new 
element in our discussion is the introduction of the lead¬ 
ing classical 0 {l/as) interactions between the dyons and 
anti-dyons as recently obtained in [52] using the classical 
“streamline” set of configurations for MM, LL pairs. We 
have assumed that the ML, LM channels are repulsive 
and opposite in sign to the streamline interaction. While 
carrying this work, this assumption has now been con¬ 
firmed numerically |35j . Another important element of 
our analysis is the one-loop measure of the dyon and anti- 
dyon moduli space, in the form proposed by Diakonov 
and Petrov |1|. It leads to a small moduli space volume 
and thus repulsive interaction at higher density, which 
however can be made much less repulsive by introducing 
correlaltions between the charges. 

On general grounds, an ensemble of instanton-dyons is 
a strongly coupled plasma, with significant correlations 
between the particles. Therefore, the statistical mechan¬ 
ics of a generic instanton-dyon ensemble is very nontrivial 


and remains unsolved. However - and this is the main ar¬ 
gument of the paper - when the plasma is dense enough 
for temperatures below Tc, it generates a large screen¬ 
ing mass M which screens the interaction. A standard 
weak coupling plasma theory, in a form similar to the 
Debye-Huckel theory is then applicable . The dimension¬ 
less 3-density of each dyon species in the regime 

considered is in the range of no/d « T^/d, in agreement 
with the qualitative arguments in Ha- 

Using it, we get a number of results concerning the de¬ 
tails of the non-perturbative gauge fields, in the temper¬ 
ature range (0.5 — 1) Tc- First , in the presence of strong 
screening the minimum of the free energy is still at the 
confining (center symmetric) value of = 1/2, with a 
vanishing Polyakov line (L) « cos{2'kv) = 0 . Second, a 
re-summation of the the linearized screening effects yields 
Debye-Huckel type corrections to the pressure and dyonic 
densities. We have also analyzed the topological suscep¬ 
tibility, the gluonic compressibility, and the electric and 
magnetic gluonic condensates in this linearized approxi¬ 
mation. 


We have calculated also the electric and magnetic 
screening masses, generated by the dyon ensemble. We 
have found that the latter are larger than the former 
in the confined phase. This is qualitatively consistent 
with the existing lattice data, which however are much 
better measured for the SU(3) gauge theory rather than 
the SU(2) one we have studied here. Finally, we have 
calculated the structure factors in the electric and mag¬ 
netic sector in the linearized screening approximation as 
well. For an estimate of the transition temperature from 

= 1/2 (confinement) to = 0 (deconfinement) we 
have switched the perturbative (GPYW) holonomy po¬ 
tential m in section pH G[ For SU(2) the transition is 
observed to take place at s/oe ~ 0.88. 

In the dyonic plasma the large spatial Wilson loops 
exhibit area law, while the spatial t' Hooft loops are 
found to be 1 modulo 0[ag) Coulomb-like self-energy 
corrections. These dual behaviors were argued in |33j 
for confining gauge theories at zero temperature. We 
found them to hold in the confining dyon-ensemble in 
the regime 0.5 < T < T^. 

Needless to say, that all these predictions can and 
should be confronted with the lattice data in the cor¬ 
responding temperature range. 

Finally, let us speculate about the dyon ensemble be¬ 
yond the validity domain of the Debye-Huckel approx¬ 
imation. First of all, strongly coupled Coulomb plas¬ 
mas are tractable by certain analytic and/or numerical 
(molecular dynamics) methods, see Refs [231 HI] for sim¬ 
ilar development. Another option is to use brut force 
numerical simulations of the dyon ensemble [36] . Quali¬ 
tatively, sufficiently strongly coupled plasmas develop ei¬ 
ther (i) correlations between particles, resembling either 
a liquid with crystal-like correlations (“molten salt”), or 
(ii) particular neutral clusters, the simplest of which can 
be the LM instantons themselves or LMLM “instanton 
molecules”. Recent (unquenched) lattice simulations in- 
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dicate that the instantons and anti-instantons recombine 
into topologically neutral molecules across the transition 
temperature [Mills]- At much higher temperature, the 
perturbative gluons dwarf all classical gauge configura¬ 
tions forcing the holonomy to zero. 

One obvious extension of this work should be into the 
large number of colors N^.. Strong correlations can ap¬ 
pear, since « l/cus ~ ^ 1. Similar mechanism, 

leading to crystallization appears to take place in dense 
holographic matter where the baryons as instantons in 
the holographic direction split into a pair of dyons and 
re-arrange in salt crystals m- 

Another obvious extension of this work is to include 


fermions, which we turn to in the second paper of the 
series [42] , 
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